Event -by- event study of DCC-like fluctuation 
in ultra-relativistic nuclear collisions 
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A method based on sliding window scheme is developed to search for patches 
in the pseudorapidity-azimuth plane, on an event-by-event basis, having unusual 
fluctuation in the neutral pion fraction which may arise due to the formation of 
Disoriented Chiral Condensates (DCC) in high energy nuclear collisions. The effi- 
ciency of the method to extract the patches and the purity of the extracted sample 
■ are studied for possible experimental situations. 
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Disoriented Chiral Condensates (DCC) have been predicted to be formed in 
high energy hadronic and nuclear collisions when the chiral symmetry is tem- 



porarily restored at high temperatures. As the matter cools and expands, the 
vacuum may relax into a state that has an orientation different from the nor- 
mal vacuum. This may lead to the formation of localized domains of DCC, 
emitting low momentum pions in a single direction in isospin space [1]. This 
would lead to event-by-event fluctuation in the number of charged particles 
and photons in a given phase space 1 , since majority of photons originate 
from 7T° decay and charged particles are mostly charged pions. It has been 
estimated that the neutral pion fraction(/), corresponding to the DCC do- 
main, follows the probability distribution P(f) = where/ = N - . 

This distribution is markedly different from the standard binomial distribu- 
tion (peaking at /= ~) for generic pion production in nuclear collisions. The 
observed Centauro(iV c / l > iV 7 ) and Ant i- Cent auro(iV c / I < iV 7 ) events in cosmic 
ray collisions [2] may be due to the formation of large domains of DCC. 



1 A section of the pseudorapidity(r/)-azimuth(</>) phase space where DCC pions 
may be localized will be called a patch. 
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There is very little theoretical guidance about the nature and probability 
of DCC formation in nucleon-nucleon and nucleus-nucleus collisions. Crude 
theoretical estimates suggest that in central lead-lead collisions at the CERN 
SPS, DCC formation probability may be <10~ 3 [3]. The results of DCC search 
in cosmic ray [4] and pp experiments [5] have been largely inconclusive. In 
nucleus-nucleus collisions, upper limits at 90% confidence level have been set 
in the analysis of data from the WA98 and NA49 experiments at the SPS [6-9] 
within the context of a simple DCC model and under various assumptions. 
For domains localized in (77, (ft) the most strict limit set so far, for the top 
5% central Pb-Pb collision events in the WA98 experiment [8], are ~10~ 2 and 
~3xlCT 3 for azimuthal domain sizes 45°-90° and 90°-135° respectively. 

Several methods have been proposed to detect the DCC formation in high 
energy hadron and heavy ion collisions and also used in the analysis of exper- 
imental data. For details see the recent review by Mohanty and Serreau [1]. 
The sensitivity of the methods has been studied using a simple DCC model 
and generating simulated data containing a mixture of generic and DCC type 
events. The Discrete Wavelet Transform method [10] is found to be sensitive 
if the fraction of DCC events is at least 2xl0~ 3 for domain size in the range 
of 40°-90° in a scenario where all the pions within the domain are of DCC 
origin. Considering only pions with px <300 MeV/c being of DCC origin, 
the sensitivity is found to be at 10~ 2 level. The fluctuation measure [11] has 
similar sensitivity, ^-measure [12] is useful if each event contains 25% pions 
from DCC source and the DCC event fraction in the sample is 10~ 2 . It is thus 
found that the present methods, including those used in the analysis of data 
from the WA98 experiment, have sensitivity in the same range as the possible 
DCC formation probability. In this letter we describe a method which allows 
direct observation of DCC-like patches by examining the event structure and 
its sensitivity can be better than 10~ 4 , being limited only by the available 
statistics. 

We use a simple DCC model where charged neutral fluctuation is introduced 
in generic events produced by the VENUS event generator [13], keeping ini- 
tially the identity of ir° intact. It is assumed that the DCC pions survive till 
freezout, and that they also remain localized in a (77, <fi) patch. To introduce 
DCC-like fluctuation the identity of charged and neutral pions is interchanged 

pairwise (ty + ty^ <-> 7r°7r°) according to the — y= probability distribution within 

2 v / 

a randomly selected (77, (p) domain. Two scenarios are considered for DCC-like 
fluctuation in the selected domain: 

• DCC-I : all the pions are of DCC origin; 

• DCC-II : only the pions having p T <250 MeV/c are of DCC origin. 

For the present study the DCC domain size is kept as one unit in pseudo- 
rapidity and 60° in azimuth. 7r°'s are allowed to decay after introducing the 



2 



a x10 3 

a 

<u 

6 10 

"8 
o 
Z 



4 




Neutral pion fraction (f) 



Fig. 1. Neutral pion fraction(/) distributions for DCC-I patches before 7r°-decay 
(solid line). The other three distributions are, after 7r°-decay, for generic events 
(dash-dot line), DCC-I patches (dashed line) and DCC-II patches (dotted line). 

DCC-like fluctuation. Considering experimental situations where one measures 
the multiplicities of charged particles (N ch ) and photons (iV 7 ) in an event, the 
neutral pion fraction may be approximately written as / = N h • ^ e con " 
sider typical nuclear systems and collision centralities such that after 7r°-decay, 
the mean multiplicities of photons and charged particles are ~400 within one 
unit of pseudorapidity and with full azimuthal coverage in the generic event. 
Such multiplicities are typical for detector systems available in SPS and RHIC 
experiments. 

Neutral pion fraction (/) distribution for DCC-I patches before 7r°-decay is 
shown in Fig. 1. The /-distribution after 7r°-decay, is also shown in Fig. 1 for 
various cases. In general the /-distributions shrink inward for both scenarios 
of DCC formation after ir° decay [14], the effect being more pronounced for 
DCC-II. It is clear that in the central part of /-distribution the study of 
DCC-like fluctuation is difficult due to overwhelming background from the 
binomial distribution of generic particle production mechanism. We therefore 
concentrate on the low-/ and high-/ regions, which will be referred to as 
charge-excess and photon-excess type fluctuations respectively. 

Different ensembles of simulated DCC-like events are generated by adding 
events having DCC-like fluctuation and generic events in different proportions, 
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reflecting the DCC formation probability. Such ensembles of events will be 
referred to as 'simulated data'. The fraction of DCC-like events (a) in various 
ensembles of simulated data varies from 2xl0~ 4 to 1. The total number of 
events in each ensemble is 10 5 . 

The search for any non-statistical fluctuation requires a reference for compari- 
son which can describe statistical fluctuation in a model-independent manner 
and is free from inherent correlations. For the analysis of experimental data 
the technique of 'mixed events' provides such a reference set where correlations 
among particles are completely destroyed. The mixed events are generated us- 
ing the standard procedure of taking large number of events (at least equal 
to the multiplicity of particles) and constructing new events having particles 
picked up randomly, one particle from one event. These events are created 
with the same multiplicity distribution as the 'simulated data', similar statis- 
tics and have similar charge-neutral correlation on a global scale. 

For the study of charged neutral fluctuation, we assume that there are detec- 
tors measuring the photon and charged particle multiplicities in overlapping 
part of (77, <fi) phase space, with or without the p T -information of particles. 
The common coverage of the two detectors is assumed to be one unit in rj and 
2n in <f>, the introduced DCC domain being fully contained within this region. 
We develop a simple but powerful 'sliding window' method (SWM) in which a 
window in azimuthal plane, say of size A<f>, is chosen in the common coverage 
of the two detectors in which the neutral pion fraction / is calculated. The 
entire azimuthal range of common coverage is scanned by continuously sliding 
the window, shifting each time by a small amount, say 8<p, to search for a 
patch having a neutral pion fraction several standard deviations away from 
the mean value. This method utilizes the full advantage of azimuthal resolu- 
tion of the detectors and allows direct observation of patches having large (or 
small) /-values. The value of 5(f) depends on the azimuthal resolution of the 
two detectors. For the present work, we have chosen 5<p=2°. A preliminary 
version of the SWM is described in [15]. 

The SWM provides a set of /-values in each event. We focus on the maximum 
and minimum /-values (/ max and / m i n ) to search for photon-excess and charge- 
excess type fluctuations. The / max and / min distribution of simulated data for 
a=0.01 and 0.5 are shown in Fig. 2 for a typical window size A0=6O° along 
with those for mixed events generated from those samples. The / max and / m i n 
distributions for mixed events are almost Gaussian in shape for lower value of 
a but become skewed for large values of a. 

It is observed that the signals of non-statistical fluctuation can be distin- 
guished only in the outer regions. We estimate the mean (//) and r.m.s. devi- 
ation (a) of the mixed event distributions. A cut on the / max and / m i n distri- 
bution of simulated data is applied by fj, ± na, where n ranges from 3 to 5, 
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Fig. 2. /max(left) and /min(right) distributions for a A<fi = 60° window size for two 
values of a in DCC-I scenario as indicated. Continuous lines : simulated events, 
dashed lines : mixed events generated from corresponding simulated events. 

positive sign being applied for / max and negative sign for / min distributions. 
The events having patches where the / max or / min values are beyond the cut 
are labeled as 'extracted'. 

The sensitivity of the SWM to extract DCC-like fluctuation is decided by 
the limit of statistical background that the mixed event sample gives for any 
given extraction procedure. The patches in mixed events having very large 
(or very small) /-values, arising due to statistical fluctuation, constitute the 
background to the study of DCC-like fluctuation. The statistical fluctuation 
depends on the multiplicity of particles in the patches and hence also on the 
window size in the SWM. 

For mixed events produced from generic events (a=0) Fig. 3(a) displays the 
fraction of events, which constitutes the background, versus window size for 
/i ± 4a cut on mixed event / max and / m i n distributions. The background is 
almost independent of window size for A<p> 60° and is found to be 4xl(r 4 and 
lx 10~ 4 respectively for photon-excess and charge-excess cases. For subsequent 
investigations we have used A^=60°. 

The background due to statistical fluctuation has also been investigated using 
the mixed events generated from different ensembles of events (different a 
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Fig. 3. Fraction of events versus window size for q=0 (top) and versus a for A(/>=60 o 
(bottom) for / max distribution with /x+4cr cut (filled circles) and for / m ; n distribution 
with //-4<r cut (filled squares) for mixed events generated from different ensembles 
of events. Filled triangles in the bottom panel represent fraction of events for / max 
distribution with fi+5a cut. 

cases), keeping the window size fixed at 60°. Fig. 3(b) shows the fraction 
of events as a function of a for // ± 4a cuts on mixed event / max and / m i n 
distributions. The background is of similar magnitude as that in Fig. 3(a) 
for lower values of a. With increasing a, the background increases and then 
decreases for very high values. This apparent decrease at very high a is due 
to the distorted shape of the corresponding mixed event distributions and the 
resulting high values of a. Fig. 3(b) also displays the background for photon- 
excess cases with jj, + 5a cuts. The background values are now an order of 
magnitude smaller, being as low as 5xlCT 5 for smaller values of a. 



In the present case of simulation study, where the character of each event 
and patch is known, we can estimate the 'background' (B) which arises from 
the following sources : (a) the extracted patch is from a generic event and 
(b) the extracted patch, in the DCC-like event, is not from the azimuthal 
region where DCC was actually introduced. Extracted patches which belong 
to the DCC-like events and are found around the same azimuthal location 
are labeled as 'extracted signal' (S). In order to characterize the SWM we 
estimate the efficiency of extraction S/N, where N is the number of DCC-like 
events present in the ensemble depending on a value, and the 'signal purity 
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Fig. 4. Efficiency and purity (expressed in percentage) versus window size for sim- 
ulated data with a=0.01 and [i±4<j cut as the extraction criteria. Left panels are 
for photon-excess cases and right panels for charge-excess cases. 

of extracted patches' S/(S+B). 

The effect of window size on the results of SWM analysis has been investi- 
gated further for the case of simulated events with a=0.01. The efficiency of 
extraction and the purity of the extracted patches are plotted for the case 
of both charge-excess and photon-excess fluctuations as a function of window 
size in Fig. 4, using fi±Aa cut as the extraction criteria. It is found that the 
efficiency and purity are not much affected by varying the window size over 
a wide range around the azimuthal domain where fluctuation was originally 
introduced. The results are affected only at very small window sizes, which, as 
discussed earlier, get dominated by statistical fluctuation due to small number 
of particles. 

Fig. 5 shows the efficiency and purity as a function of a for the entire data set 
analyzed in the present work for both charge-excess and photon-excess type 
fluctuations in DCC-I and DCC-II scenarios. The extraction criteria applied 
is a cut of fi ± na on the / max and / min distributions obtained from the corre- 
sponding mixed events for each a. Results for DCC-I scenario are presented 
for n=3,4,5 for photon-excess and for n=3,4 for charge-excess cases. 

As a function of increasing n, the signal falls very slowly whereas the back- 
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Fig. 5. Efficiency and purity versus a. The left panels are for photon-excess and the 
right panels for charge-excess fluctuations. Filled circles : 3a cut for DCC-I, filled 
square : 4a cut for DCC-I, filled triangles : 5a cut for DCC-I, open stars : 4a cut for 
DCC-II with multiplicity measurement only, filled stars : 4a cut for DCC-II with 
the addition of charged particle pr measurement. 



ground falls sharply. The efficiency is almost constant for lower values of a. The 
decrease in efficiency for larger values of a is due to distortions in the mixed 
event / m i n and / max distributions and the apparently large cr-values. The effi- 
ciency for extraction of patches for both the charge-excess and photon-excess 
cases in DCC-I scenario is found to be in the range of 10-20%, lower values of 
efficiency corresponds to higher purity. Even for a as low as 3xl0 -4 , the pu- 
rity is more than 50%, indicating that the background is less than the signal. 
It is clear that if the statistics is high, one can make even tighter cuts (n >5) 
to further improve the purity with only marginal loss in efficiency. 

For the analysis of data having fluctuation in DCC-II scenario, two approaches 
are used. First we assume that only multiplicity information is available. The 
analysis proceeds in a manner identical to what was done for DCC-I case. 
The efficiency of extracting DCC-II type patches and the purity of such ex- 
tracted patches are shown in Fig. 5 using \x ± 4a cut as the extraction criteria. 
The efficiency is now quite low for both charge-excess and photon-excess type 
patches. This is understandable as the fraction of DCC pions in the patches 
is now reduced to less than half. 
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In the second approach we assume that the pr information of charged particles 
is available. Such a possibility exists in the STAR experiment at RHIC [16]. 
The neutral pion fraction / is now calculated by taking only the charged 
particles having p T <250 MeV/c. Using the new /-values both / min and / max 
distributions are generated. The efficiency of extraction and the purity of 
extracted samples, shown in Fig. 5 using the fx ± 4a cut as the extraction 
criteria, are found to increase by a factor of about 5-10 over the case when 
only multiplicity information was used. 

In summary the SWM provides a simple and elegant method for the study of 
charge-neutral fluctuation in ultra-relativistic nuclear collisions. The sensitiv- 
ity of the method, decided by the background due to statistical fluctuation, can 
be improved to a great extent by making suitable extraction criteria so that 
the background is minimized. In the present analysis, with just 10 5 events, it 
has been possible to reduce the background to 5x 1CT 5 . This is almost two or- 
ders of magnitude better than the previous methods. At higher multiplicity, as 
is expected at the LHC experiments, background due to statistical fluctuation 
may be further reduced. 

The efficiency of extracting the DCC-like patches varies from a few percent 
to about 40% (combining the photon-excess and charge-excess cases). This 
depends on the extraction criteria, DCC formation probability and the p-r 
distribution of DCC pions. The purity of the extracted patches also depends 
on those parameters and can reach 60-80% even for low values of a. For high 
values of a this reaches close to 100%. For a given a and extraction criteria the 
purity is higher for photon excess cases than for charge-excess cases as observed 
earlier [14]. In DCC-II scenario the efficiency is much lower for charge-excess 
cases than for photon-excess cases. Considering limitations due to detector 
effects [14] it is advantageous to look for photon-excess type charge-neutral 
fluctuation. 

The SWM is a general method and can be utilized not only in azimuthal 
space but also in pseudorapidity space or even in any combined phase space 
with mult i- dimensional windows and using any suitable physical observable 
which can be computed over the window. The only important point is to slide 
the window over the acceptable phase space region to compute a number of 
values of the observable under investigation and then study their distributions. 
This allows for direct observation of unusual structures in the event and draw 
conclusions which are model-independent. 

The method developed here can be directly applied to the analysis of exper- 
imental data using mixed event as the reference set without recourse to any 
DCC model. With improved statistics, which is now-a-days available in collider 
experiments, it should be possible to reduce the background to any desired 
level and improve the sensitivity of the SWM. This should also allow one to 
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employ tight cuts which further improves the signal purity to the estimated 
level to confirm or deny the occurrence of exotic phenomena. 
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